Diatomic molecules in ultracold Fermi gases - Novel composite bosons 
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We give a brief overview of recent studies of weakly bound homonuclear molecules in ultracold 
two-component Fermi gases. It is emphasized that they represent novel composite bosons, which 
exhibit features of Fermi statistics at short intermolecular distances. In particular, Pauli exclusion 
principle for identical fermionic atoms provides a strong suppression of collisional relaxation of such 
molecules into deep bound states. We then analyze heteronuclear molecules which are expected to 
be formed in mixtures of different fermionic atoms. It is found how an increase in the mass ratio for 
the constituent atoms changes the physics of collisional stability of such molecules compared to the 
case of homonuclear ones. We discuss Bose-Einstein condensation of these composite bosons and 
draw prospects for future studies. 



PACS numbers: 34.50.-s, 03.75.Ss 
I. INTRODUCTION 

Quantum statistics for bosons and fermions lies in the 
basis of our understanding of many-body physics. In the 
case of a dilute gas of bosons, one faces the phenomenon 
of Bose-Einstein condensation (BEC), a macroscopic oc- 
cupation of a single quantum state below a certain criti- 
cal temperature, predicted by S. Bose and A. Einstein in 
1924 [JQ. The idea of BEC originally comes from the 
analogy between light and matter waves. Indeed, at suf- 
ficiently low temperatures the de Broglie wavelength of 
particles exceeds the mean inter-particle separation, and 
the wave packets of particles start overlapping. Then the 
quantum statistics comes into play and it becomes favor- 
able for bosons to occupy a single quantum state. This 
state represents a macroscopic quantum object called 
Bose-Einstein condensate, and the formation of the con- 
densate manifests itself as a phase transition accompa- 
nied by a change of some of the properties of the gas. 

Since the work of Bose p| and Einstein Q , a number of 
phenomena in interacting systems have been considered 
as manifestations of BEC: superfluidity in liquid helium, 
high-T c superconductivity in some materials, condensa- 
tion of hypothetical Higgs particles, BEC of pions, and so 
on. Bose-Einstein condensation in dilute gases has been 
first observed in 1995 in pioneering experiments with 
clouds of magnetically trapped alkali atoms at JILA 
MIT and Rice The first generation of studies of 
trapped Bose-condensed gases found a variety of spectac- 
ular macroscopic quantum effects: interference between 
two condensates, collective oscillations and their damp- 
ing, formation of quantized vortices and vortex lattices, 
atom lasers, etc. (see @ for review). 

In the last few years, the field of quantum gases is 



strongly expanding in the direction of ultracold two- 
component clouds of fermionic atoms, with the initial 
goal of achieving a superfluid Bardeen-Cooper-Schrieffer 
(BCS) phase transition. This transition requires attrac- 
tive interaction between the atoms. Then, in the sim- 
plest version of the superfluid transition, at sufficiently 
low temperatures fermions belonging to different compo- 
nents and having opposite momenta on the Fermi surface 
form correlated (Cooper) pairs in the momentum space. 
This leads to the appearance of a gap in the excitation 
spectrum and to the phenomenon of superfluidity (see, 
for example, 0). In a dilute ultracold two-component 
Fermi gas, most efficient is the formation of Cooper pairs 
due to the attractive intercomponent interaction in the s- 
wave channel (negative s-wave scattering length a). How- 
ever, for ordinary values of a, the superfluid transition 
temperature is extremely low. 

For this reason, the efforts of experimental groups have 
been focused on modifying the intercomponent interac- 
tion by using Feshbach resonances (see below). In this 
case, one can switch the sign and tune the absolute value 
of a, which at resonance changes from +oo to — oo. This 
has led to remarkable developments, in particular to the 
creation of weakly bound diatomic molecules of fermionic 
atoms on the positive side of the resonance (a > 0) 
@ M- • These are the largest diatomic molecules ob- 
tained so far. Their size is of the order of a and it reaches 
hundreds of nanometers in current experiments. Accord- 
ingly, their binding energy is exceedingly small (10 -10 eV 
or less). Being composite bosons, these molecules obey 
Bose statistics, and they have been Bose-condensed in 
JILA experiments with 40 K2 0, TT3l and in 6 Li? experi- 
ments at Innsbruck [HEI, MIT HBlH, ENS G3> and 
Rice PI. 

Nevertheless, some of the interaction properties of 
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these molecules reflect Fermi statistics of the individ- 
ual atoms forming the molecule. In particular, these 
molecules are found remarkably stable with respect to 
collisional decay. Being in the highest rovibrational 
state, they do not undergo collisional relaxation to deeply 
bound states on a time scale exceeding seconds at den- 
sities of about 10 13 cm -3 , which is more than four or- 
ders of magnitude longer than the life time of similar 
molecules consisting of bosonic atoms. 

Currently, a new generation of experiments is being set 
up for studying degenerate mixtures of different fermionic 
atoms, where one expects the formation of heteronuclear 
weakly bound molecules. In this paper we analyze elastic 
interaction between these molecules and their collisional 
relaxation to deeply bound states. The emphasis is put 
on the dependence of these properties on the mass ra- 
tio for the constituent atoms and on the comparison of 
the results with those for presently studied homonuclear 
molecules. We present a physical picture of how such 
composite bosons can exhibit features of Fermi statistics 
and show that this is the reason for their collisional sta- 
bility. 



II. FESHBACH RESONANCES AND DIATOMIC 
MOLECULES IN COLD FERMI GASES 

The studies of interacting quantum gases are actively 
pursued using Feshbach resonances, in particular for tun- 
ing and increasing the intercomponent interaction in two- 
species Fermi gases. In the vicinity of the resonance, for 
the two-particle problem one has a strong coupling be- 
tween (zero-energy) continuum states of colliding atoms 
and a bound molecular state of another hypcrfine man- 
ifold. The effective coupling strength depends on the 
detuning from the resonance (the energy difference S 
between the bound molecular state and the continuum 
states), which can be varied by changing the magnetic 
field. Thus, the scattering length becomes field depen- 
dent (see Fig. 1). 

At resonance the scattering length changes from +00 
to —00, and in the vicinity of the resonance one has the 
inequality kp\a\ ^ 1 where /cf is the Fermi momentum. 
The gas then is in the strongly interacting regime. It is 
still dilute and the mean interparticle separation greatly 
exceeds the characteristic radius of interparticle interac- 
tion R e . However, the amplitude of binary interactions 
(scattering length) is larger than the mean separation be- 
tween particles, and the ordinary mean field approach is 
no longer valid. 

For a large detuning from resonance the gas is still 
in the weakly interacting regime, i.e. the inequality 
n\a\ 3 <C 1 is satisfied, where n is the density of the gas. 
On the negative side of the resonance (a < 0), at suf- 
ficiently low temperatures of a two-component atomic 
Fermi gas one expects the BCS pairing between distin- 
guishable fermions, well described in literature |7(. On 
the positive side (a > 0) two fermions belonging to differ- 
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FIG. 1: The dependence of the scattering length on the mag- 
netic field. The symbols I, II, and III stand for the regime 
of a weakly interacting atomic Fermi gas, strongly correlated 
regime, and the regime of weakly bound molecules. At low 
enough temperatures region I corresponds to the BCS super- 
fluid pairing, and region III to Bose-Einstein condensation of 
molecules. 



ent components form weakly bound molecules. The col- 
lisional stability and Bose-Einstein condensation of these 
molecules is a subtle question, and it is the main topic of 
this paper. 

The idea of resonance coupling through a Feshbach 
resonance for achieving a superfluid phase transition in 
ultracold two-component Fermi gases has been proposed 
in Refs. [H HJ. The crossover from the BCS to BEC 
behavior attracts now a great deal of interest, in par- 
ticular with respect to the nature of superfluid pairing, 
transition temperature, and elementary excitations. It 
is worth noting that this type of crossover has been dis- 
cussed in literature in the context of superconductivity 
[2ll I22I I2H , |24| and in relation to superfluidity in two- 
dimensional films of 3 He [25t |26| . 

The description of a many-body system near a Fesh- 
bach resonance requires a detailed knowledge of the 2- 
body problem. This is a two-channel problem which 
can be described in terms of the Breit-Wigner scattering 
[27], |2£| , the open channel being the states of colliding 
atoms and the closed channel the bound molecular state 
of the other hyperfine domain. Various aspects of this 
type of problems have been discussed by Feshbach j2^ 
and Fano [3(J • In cold atom physics the idea of Feshbach 
resonances was introduced in Ref. |3l| . and optically in- 
duced resonances have been discussed in Refs. [32ll3a |. 

The two-body physics is the most transparent if one 
can omit the (small) background scattering length. Then 
for low collision energies e the scattering amplitude is 
given by [2^| : 



F(s) = - 



+ S + ijy/e' 



(1) 



where the quantity ^7 / 'y/2fi = W characterizes the cou- 
pling between the two hyperfine domains and /1 is the 
reduced mass of the two atoms. The scattering length is 
a = —F(0). In Eq. the detuning S is positive if the 
bound molecular state is below the continuum of colliding 
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atoms. Then for S > the scattering length is positive, 
and for S < it is negative. Introducing a characteristic 
length 



R* = h 2 /2fiW 



(2) 



and expressing the scattering amplitude through the rel- 
ative momentum of particles k = y/2^,e/h, Eq. takes 
the form: 



F(k) = - — 

n J- 



1 



R*k 2 +ik 



(3) 



The validity of Eq. does not require the inequality 
kR* C 1. At the same time, this equation formally co- 
incides with the amplitude of scattering of slow particles 
by a potential with the same scattering length a and an 
effective range R = —2R*, obtained under the condition 
k\R\ < 1. 

The length R* is an intrinsic parameter of the Feshbach 
resonance problem. It characterizes the width of the res- 
onance. From Eqs. and J2J we see that large W and, 
hence, small R* correspond to a wide resonance, whereas 
small W and large R* lead to a narrow resonance. The 
issue of wide and narrow resonances is now actively dis- 
cussed in literature [MIMIIIIIISIMIIIIIII 

In fact, the use of the terms "wide" and "narrow" de- 
pends on the problem under consideration. For example, 
in the unitarity limit where a — > ±oo, Eq. © shows that 
the length R* drops out of the problem under the con- 
dition kR* <C 1. In a quantum degenerate Fermi gas 
the characteristic momentum of particles is the Fermi 
momentum kp — (S^n) 1 ^ 3 . Thus, for a given R* the 
condition of the wide resonance depends on the gas den- 
sity n and takes the form k F R* < 1 [H HUMES E3 • 
Then the scattering properties of the gas are the same as 
in the case of potential scattering at an infinite scattering 
length. 

For a > and na 3 -C 1, where one has a weakly inter- 
acting gas of diatomic bosonic molecules I42j . the crite- 
rion of the wide resonance is different [38l |43| . The most 
important limitation is related to the binding energy and 
the wavefunction of the molecules. The energy of the 
weakly bound molecular state (it is certainly assumed 
that the characteristic radius of interaction R e <C a) is 
determined by thepole of the scattering amplitude 
One then finds [38LI43] that this state exists only for a > 
and under the condition 



R* < a. 
the binding energy is given by 

Eq = h 2 /2/ia 2 . 



(4) 



(5) 



Then the wavefunction of the weakly bound molecular 
state has only a small admixture of the closed channel, 
the size of the molecule is ~ a, and atom-molecule and 
molecule-molecule interactions are determined by a sin- 
gle parameter - the atom-atom scattering length a. In 



this sense, the problem becomes universal. It is equiva- 
lent to the interaction problem for the two-body poten- 
tial which is characterized by a large positive scattering 
length a and has a potential well with a weakly bound 
molecular state. The picture remains the same when the 
background scattering length ab g can not be neglected, 
although the condition of a wide resonance can be some- 
what modified |44j . 

Experimental studies are currently dealing with Fermi 
gases consisting of atoms in two different internal (hyper- 
fine) states and, accordingly, having equal masses. Most 
experiments use wide Feshbach resonances. For exam- 
ple, weakly bound molecules ^Li? and 40 Kq have been 
produced in experiments d U Efl [H E M, III 13 HI 
Il7l fl8| by using Feshbach resonances with a length R* of 
the order of or smaller than 20A. Then, for the achieved 
values of the scattering length a from 500 to 2000A, the 
ratio R* /a is smaller than 0.1. The only exception is the 
experiment at Rice with 6 Li near a narrow Feshbach res- 
onance at 543 G m| . For this resonance the length R* is 
very large and at obtained values of a the condition Q 
is not fulfilled. Below we will focus attention on the case 
of a wide Feshbach resonance. 



III. ELASTIC INTERACTION BETWEEN 
WEAKLY BOUND MOLECULES 

At temperatures sufficiently lower than the molecular 
binding energy £o and for equal concentrations of the two 
atomic components, practically all atoms are converted 
into molecules if the gas density satisfies the inequality 
na 3 -C 1 0| . This will be definitely the case at temper- 
atures well below the Fermi energy (the least of the two 
Fermi energies in the case of atoms with different masses). 
Thus, one has a molecular Bose gas and the first key 
question is the elastic interaction (scattering) between 
the molecules. It is important for evaporative cooling 
of the molecular gas to the regime of Bose-Einstein con- 
densation and for the stability of the condensate. In the 
latter respect, the elastic interaction should be repulsive, 
otherwise the Bose-condensed molecular gas undergoes a 
collapse. 

The exact solution for the molecule-molecule elastic 
scattering in the case of equal atom masses was found 
in Ref . |43| , and a detailed analysis of this problem was 
given in Ref. 0] . This was done assuming that the (pos- 
itive) atom-atom scattering length a for the interspecies 
interaction greatly exceeds the characteristic radius of 
interatomic potential: 



a > R e 



(0) 



Then, as in the case of the 3-body problem with fermions 
[3 HI 01 l50j . the amplitude of elastic interaction is 
determined only by a and can be found in the zero-range 
approximation for the interatomic potential. 

This approach was introduced in the two-body physics 
by Bethe and Peierls The leading idea is to solve 
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the equation for the free relative motion of two particles 
placing a boundary condition on the wavefunction ip at 
a vanishing interparticle distance r: 



NO' 

rip 



1 



0. 



which can also be rewritten as 



%l> oc (l/r — 1/a), r — > 0. 



(7) 



(8) 



One then gets a correct expression for the wavefunction 
at distances r ^> R e . For the case where a » R e , 
Eq. iJSJ gives a correct result for the wavefunction of 
weakly bound and continuum states even at distances 
much smaller than a. 

In this section we generalize the results on the 
molecule- molecule scattering obtained in Refs. *° 
the case of molecules consisting of fermionic atoms with 
different masses M and m (m < M). This is relevant 
for weakly bound heteronuclear molecules, such as 6 Li- 
40 K, 6 Li- 87 Sr, 40 K- 173 Yb etc., which can be formed by 
sweeping across a Feshbach resonance for the interspecies 
interaction. Elastic interaction between such molecules 
is important for understanding the physics of their Bose- 
Einstein condensation and for studying the BCS-BEC 
cross-over in Fermi mixtures. 

The ultracold limit for the molecule-molecule scatter- 
ing is realized under the condition k « 1, where k is the 
relative momentum. This is because a is approximately 
the range of interaction between the molecules. In this 
case the scattering is dominated by the contribution of 
the s-wave channel. The inequality ka <C 1 is equivalent 
to the collision energy much smaller than the molecular 
binding energy eq. Hence, the s-wave molecule-molecule 
elastic scattering can be analyzed putting the total en- 
ergy equal to — 2eq = —h 2 /^ia 2 . In the zero-range ap- 
proximation one should solve the four-body free-particle 
Schrodinger equation: 



R 



2/a 2 ] * = 0, 



(9) 



placing the Bethe-Peierls boundary condition at van- 
ishing distances between heavy and light fermions. In 
Eq. |(?JJ the distance between two given heavy and light 
fermions is ri, and r 2 is the distance between the other 
two (see Fig. 2). The distance between the centers of 
mass of these pairs is /3R, and r± = a±ri + a T r 2 ± /3R 
are the separations between heavy and light fermions in 
the other two possible heavy-light pairs. Here we intro- 
duced constants P = y/2a + a-, a+ — (i/M, a_ = fj,/m, 
and the reduced mass is fi — mM/(m + M). 

The wavefunction \1/ is symmetric with respect to 
the permutation of composite bosons and antisymmet- 
ric with respect to permutations of identical fermions: 

*(!-!, r 2 ,R) = *(r 2 ,r!,-R) 

= -^(r^r^i/^r! - r 2 ) =p (o+ - a_)R). (10) 

The Bethe-Peierls boundary condition should be placed 
for a vanishing distance in any pair of heavy and light 
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FIG. 2: Set of coordinates for the four-body problem. 



fermions, i.e. for ri — * 0, r 2 — > 0, and r± — > 0. Due to 
the symmetry it is necessary to require a proper behavior 
of ^ only at one of these boundaries. For ri — > the 
boundary condition reads: 



tt(n,r 2 ,R) -» /(r 2 ,R)(l/47rri - 1/4tto). 



(11) 



The function /(r 2 , R) contains the information about the 
second pair of particles when the first two are on top of 
each other. 

For large R the wavefunction is given by 



* « <Mri)<Mr 2 )(l - add/ PR); R > a, 



(12) 



where add is the molecule-molecule scattering length, and 
the wavefunction of the weakly bound molecule is given 
by 



4>a{r) = (rV27ra) 1 exp(— r/a). 



(13) 



Combining Eqs. (|ll|l and ljl2|l we obtain the asymptotic 
expression for / at large distances R: 

/(r 2 , R) « (2/r 2 a) exp (-r 2 /a)(l - a dd /PR); R > a. 

(14) 

For the s-wave scattering the function / depends only 
on three variables: the absolute values of r 2 and R, and 
the angle between them. An equation for the function 
/ in the case m = M is derived in Refs. |42l l43l |. The 
generalization of this approach to the case of different 
masses is straightforward and we obtain the same integral 
equation: 

j {G(|5i-Si|)[/(r',R')-/(r,R)]+ [G{\S X -S 2 \) 



R' 



(V2-l)/(r,R) 



47TO 



(15) 



The 9-dimensional vectors Si, S2, and 5i are again given 
by Si = {0,r',R}, S 2 = {r',0,-R'}, §1 = {0,r,R}, 
and G(x) = (2Tr)- 9 ^(xa/V2)- 7 / 2 K 7/2 (V2x/a) is the 
Green function of Eq. ©, with K 7 / 2 being the decay- 
ing Bessel function. The effect of different masses is con- 
tained in the expressions for the vectors S±, which now 
read ^ = {a T r'±PK', a±r'^0R', =F^r'=F(a+-a_)R'}. 
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FIG. 3: The ratio add /a, versus M/m. 



It is more convenient to make calculations in the mo- 
mentum space, transforming Eq. (|15f) into an equation 
for the function f(k, p) = J d 3 rd 3 Rf(r, R) exp(ik • r/a+ 
i(3p ■ R/a): 



E 



/(k±a T (p'-p),p')dV 



„ 2 +/3V 2 +(k ± a T (p'- p)) 2 + (k ± a±(p'+ p)) 
/(k', -p) d 3 fc' 2tt 2 (1 + k 2 + /?V)/(k, p) 



2 + fc' : 



-f3 2 P 2 



y/2 + k 2 + I3 2 p 2 + 1 



(16) 



By making the substitution /(k, p) = (<5(p) + 
,g(k, p)/p 2 )/(l + k 2 ) we reduce Eq. (|Po)l to an inhomoge- 
neous equation for the function g(k, p). For p — > this 
function tends to a finite value independent of k. The 
molecule-molecule scattering length is given by add — 
— 2ir 2 a linip^o g(k, p), and we have calculated this quan- 
tity numerically from Eq. (|16[) . 

In Fig. 3 we display the ratio add jo, versus the mass 
ratio M/m. For m — M we recover the result of Refs. 
|42l |43|: add = 0.6a. Compared to earlier studies which 
were assuming add — 2a j24[, this exact result gives al- 
most twice as small a sound velocity of the molecular 
condensate and a rate of elastic collisions smaller by an 
order of magnitude. An approximate diagrammatic ap- 
proach has been developed in Ref . [S^l , and it has led to 
add = 0.75a. 

The universal dependence of add/ a, on the mass ra- 
tio, presented in Fig. 3, can be established in the zero- 
range approximation only if M/m is smaller than approx- 
imately 13.6. Calculations then show the absence of four- 
body weakly bound states, and for M/m ~ 1 the behav- 
ior of / suggests a soft-core repulsion between molecules, 
with a range ~ a. For the mass ratio larger than the crit- 
ical value, the description of the molecule-molecule scat- 
tering requires a three-body parameter coming from the 
short-range behavior of the three-body subsystem con- 
sisting of one light and two heavy fermions [471 l4q . 



IV. BORN-OPPENHEIMER APPROACH 

In the case of large M/m one can develop a transpar- 
ent physical picture of the molecule-molecule interaction 
by using the Born-Oppenheimer approximation |53j . In 
this case heavy fermionic atoms are moving slowly in the 
field produced by the exchange of fast light atoms. The 
adiabatic behavior assumes that the four-body problem 
can be split into two parts. First, we calculate the wave- 
functions and binding energies of two light fermions in 
the field of two heavy atoms fixed at their positions Xi 
and X2. As the motion of heavy atoms is slow, the bound 
states of the light atoms adiabatically adjust themselves 
to a given value of the separation between the heavy ones, 
x = |xi — X2I. The sum of the corresponding binding en- 
ergies gives the effective interaction potential U e f / (x) for 
the heavy fermions. An obvious second step is to solve 
the Shrodinger equation for their motion in this effective 
potential. 

For x > a, there are two bound states of a light atom 
interacting with a pair of fixed heavy atoms. The wave- 
functions of these states are given by 

V>X, ± (V) OC € — : ± : : , (17) 



|y- x i| 



|y- x 2i 



where y is the coordinate of the light atom. The func- 
tions k±(x) > are determined by the requirement that 
for y — * Xi ; 2 the wavefunction (|17fl satisfies the Bcthc- 
Peierls boundary condition ((SJ in which r is substituted 
by |y — Xi^l- This gives algebraic equations: 

n ± (x) T e~ K±{x) = x/a. (18) 

The energies of the bound states are given by 

e ± (x) = -h 2 K 2 ± (x)/2mx 2 . (19) 

Since the light fermions are identical, their two-body 
wavefunction can be constructed as an antisymmetrizcd 
product of ip* t + and ip x ,-'. 

^x(yi,y2) oc V'x,+(yi)V'x,-(y2) - ^x,+(y2)^ x ,-(yi), 

(20) 

and their total energy is U e ff(x) = £+(x) + e_(x). Note 
that the wavefunction (|20|l is antisymmetric with respect 
to the transformation x — » x. This means that solving 
the Shrodinger equation for the heavy fermions, 



--V 2 x + U eff (x) 



(21) 



we have to look for a solution that is symmetric with 
respect to their permutation (in our case it corresponds 
to the s-wave scattering). Then the total four-body wave 
function *(xi, x 2 , yi, y 2 ) oc x(x)^ x (yi, y 2 ) has a proper 
symmetry. 

Solving Eqs. (fT5f) at distances x > a gives a purely 
repulsive potential monotonically decreasing with x. At 
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distances x a it behaves as 

U e ff(x) ~ —2eq + (h 2 /2m) exp(—2x/a)/ax; x 3> a. 

(22) 

We see that the ratio of the a;-dependent part of the ef- 
fective repulsion to the kinetic energy term in Eq. 121|) 
increases with M/m. This is consistent with exact cal- 
culations based on Eqs. I|15|l and (|16[1 . 

The Born-Oppenheimer approach fails at distances 
close to a as the energy £_(x) vanishes at x = a, and 
the atom in the state "— " moves slower than the heavy 
atoms. This leads to a contradiction with the adiabatic 
approximation. However, we see that the light atom in 
the state "— " becomes essentially delocalized and at dis- 
tances x <C a the problem reduces to a three-body prob- 
lem with two heavy fermions and one light fermion in 
the state "+". As has been mentioned, this problem is 
characterized by a critical mass ratio « 13.6 above which 
the behavior of the system drastically changes. For our 
four-body problem this means that the effective poten- 
tial at these distances becomes attractive and can sup- 
port bound states. For an overcritical mass ratio, we also 
expect a resonance dependence of molecule-molecule col- 
lisions on the scattering length or short-range parameters 
of the system. 

Finally, it is worth noting that for a very large M/m 
these resonances should be very narrow, since the repul- 
sive effective potential at x > a is very strong. Therefore, 
the molecule-molecule scattering length should (on aver- 
age) increase logarithmically with M/m. Obviously, the 
effective repulsive barrier is important for the analysis of 
the inelastic losses in the molecular gas. We discuss this 
question in more detail in Sec. VI. 



V. SUPPRESSED COLLISIONAL RELAXATION 

The most exciting physics with weakly bound bosonic 
molecules of fermionic atoms is related to their colli- 
sional stability. Actually, these are molecules in the high- 
est rovibrational state and they undergo relaxation into 
deeply bound states in molecule-molecule (or molecule- 
atom) collisions. The released binding energy of a deep 
state is ~ ft 2 /2p,R 2 . It is transformed into the kinetic en- 
ergy of particles in the outgoing collisional channel and 
they escape from the sample. Thus, collisional relax- 
ation determines the lifetime of the Bose gas of weakly 
bound molecules and is therefore a crucial process. Sev- 
eral experiments show that such molecules consisting of 
bosonic 87 Rb [5J, |53 and 133 Cs atoms 56] undergo a 
rapid collisional decay. On the other hand, first observa- 
tions of weakly bound molecules Li2 and K2, consisting 
of fermionic atoms 0, 0, ITTl |45| , showed that they are 
long-lived at densities ~ 10 13 cm~ 3 . 

We now arrive at the point where quantum statistics 
of composite bosons comes into play for weakly bound 
molecules of fermionic atoms. First, we discuss the case 
of M/m ~ 1, where the interaction between molecules 



is characterized by a soft-core repulsion which does not 
significantly influence their behavior at short distances. 
Clearly behaving themselves as point-like bosons at large 
intermolecular distances, these molecules "start remem- 
bering" that they consist of fermions when the inter- 
molecular separation becomes smaller than the molecule 
size (~ a). As was explained in Refs. the k- ev rea- 

son for the remarkable collisional stability of such weakly 
bound molecules of fermionic atoms is Fermi statistics for 
the atoms in combination with a large size of the molec- 
ular state (small momenta of bound fermionic atoms). 
The physical picture is the following. The binding en- 
ergy of the molecules is Eq = h 2 /2fia 2 and their size is 
close to a. The size of deep bound states is of the or- 
der of the characteristic radius of interaction R e <C a. 
Hence, the relaxation requires the presence of at least 
three fermions at distances ~ R e from each other. As two 
of them are necessarily identical, due to the Pauli exclu- 
sion principle the relaxation probability acquires a small 
factor proportional to a power of (qR e ), where q ~ 1/a 
is a characteristic momentum of the atoms in the weakly 
bound molecular state. 

The inequality a^$> R e allows one to obtain the depen- 
dence of the molecule-molecule wavefunctions at short 
interparticle distances on the two-body scattering length 
a and thus to establish a dependence of the relaxation 
rate on a. In Refs. (52, ^| we have done this for m = M 
and found a strong decrease of relaxation with increas- 
ing a. Here we generalize these findings to the case of 
m =/= M, assuming that the mass ratio M/m is smaller 
than the critical value 13.6 and short-range physics does 
not influence the a-dependence of the relaxation rate. 

The relaxation is essentially a three-body process 
which occurs when three fermions approach each other 
to a short distance of the order of R e . The fourth parti- 
cle does not participate in the sense that the configura- 
tion space contributing to the relaxation probability can 
be viewed as a system of three atoms at short distances 
~ R e from each other and a fourth atom separated from 
this system by a large distance ~ a. In this case and also 
for any hyperradius of the three-fermion system p <C a, 
the four-body wavefunction decomposes into a product: 

* = v ( z )^ 3 \ P ,n) 1 (23) 

where is the wavefunction of the three-fermion sys- 
tem, fi is the set of hyperangles for these fermions, z is 
the distance between their center of mass and the fourth 
atom, and the function rj(z) describes the motion of this 
atom. This function is normalized to unity and aver- 
aging over the motion of the fourth atom we see that 
the probability of relaxation is determined solely by the 
short-range behavior of the function ty^(p, Q). 

In the case of fermionic atoms with different masses one 
has two possible choices of a three-body subsystem out 
of four fermions. The most important is the relaxation in 
the system of one atom with the mass m and two heavier 
atoms with masses M . 



7 



The function for p <C a follows from the three- 
body Schrodinger equation with zero total energy and, 
hence, can be written as \E , *- 3 - ) = A(a)tp(p, f2), where the 
function ip is independent of a. Assuming that the in- 
elastic amplitude of relaxation is much smaller than the 
amplitude of elastic scattering, the dependence of the 
relaxation rate constant on a is related only to the in- 
dependence of the function ty( 3 > and, hence, is given by 

a rel cx \A(a)\ 2 . (24) 

For finding a re i(a) it is sufficient to consider distances 
where R e <§C p *C a. At these distances we have ip = 
$(fi)/9" _1 0, and the function 5' (3) becomes 

*( 3 ) = A(a)^(il)p u - 1 . (25) 

The prefactor A(a) can be found from the solution of the 
four-body problem in the zero-range approximation. 

Quite elegantly, the a-dependence of A(a) can be de- 
termined from the following scaling considerations. The 
scattering length a is the only length scale in our prob- 
lem and we can measure all distances in units of a. For 
example, rescaling p — ap' and z = az' , taking into ac- 
count that rj(z) is normalized to unity, and using \I/( 3 ) 
from Eq. i|25[l we obtain that the a-dependent coefficient 
in Eq. I|23|l equals A(a)a v ~ 5 / 2 . By applying the same 
rescaling to Eq. (fT^fl with the account of Eq. lfT^|) we 
see that the same coefficient should be proportional to 
a~ 3 . Therefore, A(a) oc a _!y_1 / 2 and a re i oc a~ s , where 
s=2v + l. 

The exponent s depends on the mass ratio and on the 
symmetry of the three-body wave function \I/( 3 ). The 
leading relaxation channel at large a corresponds to the 
smallest v. The analysis of the three-body wavefunctions 
for the s-wave and p-wave atom-molecule collisions has 
been performed in Ref. 48]. The p-wave symmetry in 
the system of one light and two heavy fermions provides 
the smallest value of v which in the interval — 1 < v < 2 
is given by the root of the function 0] 

. . i/(i/ + 2) ttv ^shi7COs(^7 + 7) — sm(uj) 

A(i-M = cot -f- ~ 

v+1 2 [y + 1) sin 7COS7sin(7rz//2) 

(26) 

where 7 = arcsin [M / (M + m)] . For equal masses we 
recover s — 2v + 1 w 2.55 obtained in Refs. 0>E1, and 
it slowly decreases with increasing the mass ratio (see 
Sec. ED. 

The absolute value of the relaxation rate constant is 
determined by the behavior of the three-body system at 
short interatomic distances. Assuming that the short- 
range physics is characterized by the length scale R e and 
the energy scale h 2 /mR 2 , we can restore the dimensions 
and write: 

a rc i = C(hR e /m)(R e /a) s , (27) 
where the coefficient C depends on a particular system. 
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FIG. 4: Two-body decay rate a Te i of a 6 Li2 molecular con- 
densate a a function of interatomic scattering length near the 
6 Li Feshbach resonance at 834 G. Solid line: least square fit, 
a re i oc a _1 ' 9±0 ' 8 . Dashed line, theory: a rei oc a -2 ' 55 . The 
theoretical relaxation rate has been normalized to the exper- 
imental value at a = 78 nm. 

Experimental studies of bosonic molecules produced in 
a Fermi gas by using wide Feshbach resonances at JILA, 
Innsbruck, MIT, and ENS are well described within the 
presented theoretical approach. The JILA [hoi [Til ll^ . 
Innsbruck 0, E H, MIT [H ENS and 
Rice 0] results show a remarkable collisional stability 
of weakly bound molecules K2 and Li2 consisting of two 
fermionic atoms. At molecular densities n ~ 10 13 cm' 3 
the lifetime of the gas ranges from tens of milliseconds to 
tens of seconds, depending on the value of the scattering 
length a. A strong decrease of the relaxation rate with 
increasing a, following from Eq. (|27|) . is consistent with 
experimental data. The potassium experiment at JILA 
|lfj| and lithium experiment at ENS [131 give the relax- 
ation rate constant add oc a~ s , with s m 2.3 for K2, and 
s w 1.9 for Li2, in agreement with our theory (s w 2.55) 
within experimental uncertainty. The experimental and 
theoretical results for lithium are shown in Fig. 4. The 
absolute value of the rate constant for a 6 Li2 conden- 
sate is add ~ 1 x I0~ 13 cm 3 /s for the scattering length 
a ~ 110 nm. For K2 it is by an order of magnitude higher 
at the same value of a which can be a consequence 
of a larger value of the characteristic radius of interaction 
R e - 



VI. COLLISIONAL RELAXATION OF 
MOLECULES OF HEAVY AND LIGHT 
FERMIONS 

The analysis in the previous section also holds for col- 
lisional relaxation of weakly bound molecules of heavy 
and light fermionic atoms (M/m 1). The most ef- 
ficient is again the relaxation process in which one light 
and two heavy fermions approach each other to distances 
~ R e , and the fourth (light) atom is separated from this 
three-body system by a large distance ~ a. 

However, for a large mass ratio the behavior of the 
three-body system of one light and two heavy atoms 
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changes, and the physical picture of the relaxation pro- 
cess is drastically modified. When the heavy atoms are 
separated from each other by a distance x <C a, the 
light atom mediates an effective 1 /x 2 attraction between 
them. This follows from the Efimoy picture of effective 
interaction in a three-body system [47| and can be illus- 
trated in the Born-Oppenheimer approximation, relying 
on the results in Sec. II VI 53]. For x <§; a, the light atom 
is in the bound state "+" with two heavy atoms, and 
from Eqs. ||T5J|, (|19|l we then obtain an effective potential 
U e ff — e + w — 0.16k 2 /mx 2 . This attraction competes 
with the Pauli principle which in terms of effective inter- 
action manifests itself in the centrifugal 1/x 2 repulsion 
between the heavy atoms. The presence of this repulsion 
is clearly seen from the fact that in the state "+" the per- 
mutation of the heavy fermions does not change the sign 
of the light-atom wavefunction -0 Xj+ (y) given by Ea. (ll7|l . 
Hence, the total wavefunction of the three-body system 
V-'x.+ (y)x( x ) is antisymmetric with respect to this per- 
mutation only if the wavefunction of the relative motion 
of heavy atoms x( x ) changes its sign. Therefore, x( x ) 
contains partial waves with odd angular momenta, and 
for the lowest angular momentum (p-wave) the centrifu- 
gal barrier is U c (x) = 2h 2 /Mx 2 . For comparable masses 
it is significantly stronger than U e ff(x). Thus, we have 
the physical picture of the previous section: the Pauli 
principle (centrifugal barrier) reduces the probability for 
the atoms to be at short distances and, as a consequence, 
the relaxation rate decreases with increasing the atom- 
atom scattering length a. 

The role of the effective attraction increases with 
M/m. As a result, the decrease of the relaxation rate 
with increasing a becomes weaker. In Fig. 5 we present 
the dependence of the exponent s in Eq. (|27|) on the mass 
ratio, obtained from Eq. (|26|l which follows from the ex- 
act solution of the three-body problem. The exponent 
s continuously decreases with increasing M/m and be- 
comes zero for M/m — 12.33. In the Born-Oppenheimer 
picture this means that at this point one has a balance 
between the mediated attraction and the centrifugal re- 
pulsion. 

A further increase in M/m makes s negative and it 
reaches the limiting value s = — 1 for the critical mass 
ratio M/m = 13.6. Thus, in the range 12.33 < M/m < 
13.6 the relaxation rate increases with a. 

For an overcritical mass ratio M/m > 13.6 we have 
a well-known phenomenon of the fall of a particle to 
the center in an attractive 1/x 2 potential J28j. As has 
been mentioned in Sec. 11111 in this case the shape of the 
wavefunction at distances of the order of R e can signif- 
icantly influence the large-scale behavior and a short- 
range three-body parameter is required to describe the 
system. 

The fact that for M/m > 12.33 the relaxation rate is 
not decreasing with an increase in the atom-atom scatter- 
ing length a, does not mean that the relaxation is not re- 
duced compared to the case of molecules of bosonic atoms 
with the same mass ratio. For a large M/m, at distances 
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FIG. 5: The dependence of the exponent s = 2u+l in Eq. I27H 
on the mass ratio M/m. 

x larger than a the molecules interact via a strong re- 
pulsive potential U e ff(x) given by Eq. I|22|l . Away from 
possible resonances in molecule-molecule collisions, this 
barrier reduces the amplitude of the wavefunction at in- 
termolecular distances ~ a and, hence, leads to a sup- 
pression of the relaxation rate. Estimating the corre- 
sponding tunneling probability P in the WKB approach, 
we obtain P oc exp(— B^JM/rn) where the coefficient B 
is of order unity. 

We thus see that one also expects the suppression of 
collisional relaxation for weakly bound molecules of light 
and heavy fermionic atoms. However, the suppression 
factor is independent of the atom-atom scattering length 
a and is governed by the mass ratio M/m. The mech- 
anism of this suppression originates from Fermi statis- 
tics for the light atoms, which leads to a strong repul- 
sion between molecules at large intermolecular distances. 
This result should have implications for experiments on 
the formation of molecules in mixtures of 6 Li with much 
heavier fermionic atoms, for example 87 Sr (M/m w 14.5) 
or 173 Yb (M/mw 29). 

VII. MOLECULAR BEC AND PROSPECTS 
FOR MANIPULATIONS WITH WEAKLY 
BOUND MOLECULES 

The suppression of the relaxation decay rate of weakly 
bound molecules of fermionic atoms has a crucial con- 
sequence for the physics of these molecules. At real- 
istic temperatures the relaxation rate constant add is 
much smaller than the rate constant of elastic collisions 
& r Ka dd VT, where vt is the thermal velocity. For exam- 
ple, for the Li2 weakly bound molecules at a tempera- 
ture T ~ 3/iK and a ~ 800A, the corresponding ratio 
is of the order of 1CP 4 or 10~ 5 . This opens wide pos- 
sibilities for reaching BEC of the molecules and cooling 
the Bose-condensed gas to temperatures of the order of 
its chemical potential. Long-lived BEC of weakly bound 
molecules has been recently observed for 40 K2 at JILA 



9 



[Illii and for 6 Li 2 at Innsbruck [HG3, MIT [II EH, 

ENS |17|| , and Rice 13 • Measurements of the molecule- 
molecule scattering length confirm the result add = 0.6a 
with accuracy up to 30% [Til Il7j . This result is also 
confirmed by recent Monte Carlo calculations [52J of the 
ground state energy in the molecular BEC regime. 

Remarkable achievements in the physics of ultracold 
Fermi gases attract a lot of attention and draw fascinat- 
ing prospects for future studies. The prospects are to 
a large extent related to a very long lifetime of weakly 
bound bosonic molecules of fermionic atoms, which al- 
lows interesting manipulations with these molecules. 
Arranging a deep evaporative cooling of their Bosc- 
condensed gas to temperatures of the order of the chem- 
ical potential, one can then convert the molecular BEC 
into fermionic atoms by adiabatically changing the scat- 
tering length to negative values. This provides an addi- 
tional cooling, and the obtained atomic Fermi gas will 
have extremely low temperatures T ~ 10 _2 Tp which 
can be below the BCS transition temperature j5g. At 
these temperatures one has a very strong Pauli blocking 
of elastic collisions and expects the collisionless regime 
for the thermal cloud, which is promising for identifying 
the BCS-paired state through the observation of collec- 
tive oscillations or free expansion 0, 0, [&j| [§3, E3 • 

Recent experiments at Innsbruck, ENS, JILA and MIT 
have explored the BEC-BCS crossover regime and shown 
a remarkably reversible behavior of the gas from the BEC 
region to the BCS re gion and back over duration of sec- 
onds 0, IT3I ITfiL Il71 l62l|. Experiments have also mea- 
sured the pai ring gap in the crossover region using RF 
spectroscopy |63(. The energy gap varies between the 
binding energy of the molecules in the region a > 0, and 
a density dependent (and much smaller) value in the re- 
gion a < where no stable two-body bound state ex- 
ists. Further experiments at JILA and MIT have proven 
that fermionic atom pairs present in the strongly corre- 
lated regime at a < behave as Bose-condensed pairs 
[l2l Hil l62| . For revealing the nature of this superfluid 



pairing it is promising to make experiments with rotat- 
ing Fermi gases in the strongly interacting regime, 
where one expects the formation of vortices and vortex 
lattices [65l |6(|. Another idea is to look directly at a 
signature of the long range order in the superfluid phase 
by an interference experiment as suggested in [67l ]. 

It will also be interesting to transfer the weakly 
bound molecules of fermionic atoms to their ground ro- 
vibrational state b y u sing two-photon spectroscopy, as 
proposed in Ref. [6£| for molecules of bosonic atoms. 
Long lifetime of weakly bound molecules of fermionic 
atoms at densities ~ 10 13 cm~ 3 should provide a much 
more efficient production of ground state molecules com- 
pared to the case of molecules of bosonic atoms, where 
one has severe limitations on achievable densities and 
lifetimes. One can then extensively study the physics 
of molecular Bose-Einstein condensation. Moreover, het- 
eronuclear molecules which are supposed to be formed in 
mixtures of different fermionic atoms, can be polarized 
by an electric field. One then gets a gas of dipoles inter- 
acting via anisotropic long-range forces, which drastically 
changes the physics of Bose-Einstein condensation 69] . 
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